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Abstract 



We discuss the evaluation of power corrections to hard scattering and decay pro- 
cesses for which an operator product expansion is applicable. The Wilson coefficient 
of the leading-twist operator is the difference of two perturbative series, each of which 
has a renormalon ambiguity of the same order as the power corrections themselves, but 
which cancel in the difference. We stress the necessity of calculating this coefficient 
function to sufficiently high orders in perturbation theory so as to make the uncertainty 
of the same order or smaller than the relevant power corrections. We investigate in 
some simple examples whether this can be achieved. Our conclusion is that in most 
of the theoretical calculations which include power corrections, the uncertainties are 
at least comparable to the power corrections themselves, and that it will be a very 
difficult task to improve the situation. 
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1 Introduction 



In this paper we address the problem of controlling power corrections in effective theories. As 
an example consider e + e~ annihilation into hadrons, for which the cross section is described 
by a perturbation series, computed at the parton level, plus power corrections which are 
proportional to the condensates of higher dimensional operators. In practice the value of 
the gluon condensate is obtained by comparing the experimental value of some quantity 
derived from R e + e -(Q 2 ) to its theoretical expression. This parameter is then used to predict 
many other physical quantities, such as form factors and decay constants. Given that only 
a few (typically one or two) terms of the perturbative series are known, and that the series 
are plagued by renormalon ambiguities, which are of the same order as the contribution 
from the condensate, one may wonder whether the value of the condensate is really known 
to sufficient accuracy to be used in other processes where power corrections are important 
for the theoretical predictions. This problem is not limited to the gluon condensate, but 
is also present for other important parameters of effective theories, such as the binding 
energy (A) and kinetic energy (Ai) of the Heavy Quark Effective Theory (HQET), and the 
matrix elements of higher-twist operators in deep inelastic scattering (DIS), such as those 
for the Gross-Llewellyn Smith and Bjorken sum rules. We argue that the problem is not 
solved at present, and that the uncertainties in the determination of these non-perturbative 
parameters are seriously underestimated. 

The computation of power corrections requires the evaluation of the matrix elements of 
higher-twist or higher-dimensional operators. In addition, however, it also requires the cal- 
culation of the Wilson coefficient functions to sufficiently high order of perturbation theory 
for the cancellation of "renormalon ambiguities" to be under control Q. The reason for this 
requirement is that these ambiguities are of the same order as the power corrections. Since, 
in the calculations performed up to now, only the first few terms of the perturbation series 
are known, it is not possible to check that the remaining terms are indeed negligible. By 
studying some simple examples, we will show that the knowledge of only a few terms is, in 
general, insufficient to control the power corrections. Although the results rely on some ap- 
proximations which we are forced to adopt in these examples, it is likely that our conclusions 
will remain valid in general, and that in most further theoretical effort is needed. 

Our assumption throughout this paper is that one is attempting to evaluate the first power 
corrections (i.e. the next-to- leading twist contributions 0) with an uncertainty that is smaller 

1 We will show in section [2] that, in predictions for physical quantities, the Wilson coefficient functions 
can be written as linear combinations of two (or more) perturbation series, each of which has a high order 
behaviour such that its Borel transform has a renormalon singularity. The residue of the singularity cancels 
in the combination, however. 

2 We will frequently misuse the expression "higher-twist" to mean generic power corrections, and not just 
those to light-cone dominated quantities, such as deep inelastic structure functions. 
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than these corrections. The discussion can be readily generalized to include higher-order 
power corrections, with an increase in technical complexity, but following the same concep- 
tual principles. 

Renormalon singularities, and their implications for theoretical predictions, have been stud- 
ied for some time now, and many of the ingredients of our discussion below can be found 
in refs.|]], and in the work of Mueller [[| in particular. Recently the effects of renormalons 
in predictions for the spectroscopy and decays of heavy quarks, obtained using the heavy 
quark expansion, have been studied 0-0. We summarize in section 2 the picture which has 
emerged from these references for the appearance and cancellation of renormalon ambigui- 
ties, without reproducing the underlying arguments and derivations []. Our aim in this paper 
is to examine critically the procedure necessary to calculate the power corrections, and in 
particular to investigate whether higher-twist effects are sufficiently under control that their 
inclusion reduces the uncertainties in the theoretical predictions. Our conclusion is that to 
reach a precision of the order of the power corrections is likely to be a formidable task, and 
that the main limitation comes from the truncation of the perturbation series for the Wilson 
coefficients at low orders. The presentation below extends and clarifies that of our earlier 
paper g|. 

We particularly wish to stress that, although some of the examples presented below are given 
using the lattice spacing as the cut-off, the problems discussed in this paper are completely 
general, and are not due to some peculiarity of the lattice regularization. An example, 
where the same problems are encountered as in the lattice theory, is given by the zero recoil 
inclusive sum rules || [TO]- In this approach, in order to derive a bound on Ai, an ultraviolet 
cut-off, A, is introduced on each side of the relation between the time-ordered product of 
two currents saturated with hadronic states and the corresponding quantity computed on 
quark states using the HQET. It is not surprising that perturbative corrections ~ a s (A) A 2 
appear in the bound for Ai, essentially eliminating the predictive power of the approach (1^ . 



The appearance of such power divergences, and the consequent loss of precision when they 
are subtracted in a low order of perturbation theory, are general features in the evaluation 



of power corrections ||. A more detailed discussion of this point will be given in subsec. [272 . 

The plan of the remainder of the paper is as follows. In the next section we briefly review 
the appearance of renormalon ambiguities in the Wilson coefficient functions of operator 
product expansions. The matrix elements of the higher twist or higher dimension operators 
have to be evaluated non-perturbatively. This can be done by comparing the theoretical 
expression for a physical quantity to its experimental value (where this is known) or by 
some non-perturbative method (such as lattice simulations). We show that in both cases 
the evaluation of power corrections requires the calculation of the perturbation series to 



3 It should be stressed that although these arguments, based on the renormalization group, analyticity, 
and/or explicit calculations in the large Nf limit (Nf is the number of light quark flavours), are compelling, 
they nevertheless do not constitute a formal proof. 
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sufficiently high orders for the cancellation of renormalon ambiguities to be under control 
(sections |2J] and [T2] respectively). In sections |3|-|5] we study some simple examples, in 
order to investigate numerically the precision that might be reached in evaluating power 
corrections. We start by considering a toy model, which contains many of the general 
features concerning the appearance and cancellation of renormalon singularities in operator 
product expansions (section |3|); we then proceed to the mass of a heavy quark, or equivalently 
the binding energy A (section f|), and to the determination of the gluon condensate and its 
use in phenomenological applications (section [5]). Finally section ^contains our conclusions. 

2 Power Corrections 

In this section the cancellation of renormalon ambiguities in the evaluation of hard scattering 
and decay processes is discussed [[[J. 

Consider an operator P whose matrix element (f\P\i) contains the non-perturbative effects 
for some physical process V. In general P is non-local, for example it may be the T-product 
of two electromagnetic or weak currents at small separations (as in the e + e~ annihilation 
cross-section or weak decays), or almost light-like separations (as in deep inelastic structure 
functions). In these cases P is expanded as a series of local operators, whose coefficients 
decrease as powers of the separation. In some important applications to heavy quark physics, 
QCD composite operators (represented by P) containing the field of the heavy quark are 
expanded in terms of local operators of the HQET, with coefficients that decrease as inverse 
powers of the mass of the heavy quark. 

The expansion of P (or the Fourier transform of P) in terms of local operators Oi(/i), (^(/u) 
etc., renormalized at a scale /i, takes the form: 

vmq 2 ) = ( f l P(Q 2 ) | i ) = dCQVAi 3 ) ( / 1 Oi(a0 I i )+ C2iQ Q^ 2) ( / 1 om | i )+o 

(1) 

where n > 1 and the coefficient functions Cj are independent of the states \i) and |/). Q is a 
large momentum scale; for example it may be the centre of mass energy in e + e~ annihilation, 
the momentum transfer in deep inelastic scattering or the renormalization scale (of the order 
of the mass of the heavy quark) of the local operator P in heavy quark matrix elements. 
For clarity of notation, throughout this paper we suppress the dependence of the coefficient 
functions on the coupling constant, a s (Q 2 ). We assume here that there is only one operator 
in each of the first two orders of the expansion. If this is not the case, then there is an 
additional mixing of operators, which requires only a minor modification of the discussion 
below. We will therefore not consider this possibility further. The final term of 0(1/Q n+P ) 
in eq.(|ID represents the contributions of operators of even higher dimension which will not 
be discussed here. 
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In some cases the operator O2 is protected from mixing under renormalization with Oi, 
because of the presence of some symmetry. An important example of this in heavy quark 
physics is the chromomagnetic operator ha^G^h (where h is the field of a static quark 
and G^ v is the gluon field strength tensor), which cannot mix with the lower dimensional 
operator hh because it has a different spin structure. In such cases the corresponding problem 
of renormalon singularities does not arise. For the remainder of this paper we assume that 
O2 and 0\ have the same quantum numbers, so that they can mix under renormalization. 
Another important exception to our general discussion is the difference of matrix elements of 
the kinetic energy operator taken between different hadronic states. In this case the higher 
dimensional operator, hD 2 h, can mix with the lower dimensional one, hh, but as the latter 
is a conserved current it has the same matrix element between all single hadron states. Thus 
the corresponding renormalon ambiguity cancels in the difference of matrix elements. In 
all of these exceptions the matrix elements of the higher dimensional operators (or linear 
combinations of matrix elements) are also the leading contribution to some physical quantity 
(such as the B*-B mass splitting in the case of the chromomagnetic operator). 

The usefulness of the operator product expansion in eq.(fj]) comes from the fact that the 
non-perturbative effects in the process V are contained in the matrix elements of the local 
operators Oj, whereas the coefficient functions Cj are calculable in perturbation theory. This 
fundamental property is threatened, however, by the presence of renormalon singularities in 
the Borel transform of the perturbation series of the coefficient functions, as we now explain. 

We start by assuming that the operators on the right-hand side of eq.([TJ) are renormalized in 
some scheme based on the dimensional regularization of ultraviolet divergences, such as the 
MS scheme. Then the perturbation series for the coefficient function C% is divergent, and 
moreover is not Borel-summable, due to the presence of (infra-red) renormalon singularities 
in its Borel transform. Hence the "sum" of this series is not unique, the ambiguity being 
of 0(1/Q n ). This ambiguity is cancelled by that in the matrix element of the operator O2, 
which arises as a result of an ultraviolet renormalon singularity in the Borel transform of the 
perturbation series of its matrix elements. This implies that renormalization schemes based 
on dimensional regularization do not define higher-twist or higher-dimensional operators, 
such as 2 , unambiguously, and alternative definitions have to be used. We now consider 
the two possible alternative approaches in turn. In the first of these, the problem of the 
ambiguity in the matrix elements of O2 is eliminated by directly relating two, or more, 
physical quantities to which the matrix element of O2 contributes (see section |2.1|) . In the 
second, the operator O2 is defined using a hard ultraviolet cut-off and its matrix element is 
evaluated using some non-perturbative method, such as lattice simulations (see section |2.2|) . 
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2.1 Defining The Matrix Elements of O2 Directly in Terms of 
Physical Quantities 



In this subsection we show how to eliminate the renormalon ambiguity from the matrix 
elements of the higher-twist operator O2 by sacrificing the possibility of making a theoretical 
prediction for one process, Vfi(Q 2 ) say. This is achieved by using its experimentally measured 
value, which we denote by Vff p (Q 2 ). It is assumed here that the matrix element of 0\ 
is already known from some non-perturbative calculation, or from a second experimental 
measurement, or because, as happens in some important cases, it is a conserved operator 
or even the identity operator. We now show that such a procedure requires the control of 
the difference of two series which are not Borel-summable up to a precision of 0(1/Q n ). 
Imagine that we have computed the series for C\ and C2 in perturbation theory up to, and 
including, the term of 0(a k (Q)). We take k such that the coefficients of the series for C\ are 
not already diverging because of the renormalon singularity. Then we can define the matrix 
element of O2 by the relation 

v}7(Q 2 ) = c[ k \Q 2 /fi 2 ) ( / 1 om I i > + C ' ){Q ^ 2) ( f I of\^) 1 i ) , (2) 

where the superscript (A;) on the coefficients denotes the fact that the perturbative series 
for C\ and C2 have been truncated at 0(a*(Q)), and in the case of the matrix element 
of the higher-twist operator (/|02 that its value, derived from a physical measure- 

ment combined with a perturbative calculation, depends on the order k. The definition of 
{f\02 k \n)\i} } given in eq.(|2|) is, of course, totally unambiguous, although it does depend on 
Q and on the choice of process Vfi(Q 2 ), as well as on k. Different choices of these parameters 
can change the value of the matrix element of O2 by terms of 0(Aq CD ) (or even larger terms 
if k is too small), i.e. by an amount which is of the same order as the expected value of the 
matrix element itself. We assume throughout this paper that the value of a s is known to 
the required accuracy 0. All the elements in eq.(§) are known except for the matrix element 
of O2: Vjl P {Q 2 ) from experimental measurement, {f\0\{ji)\i) by symmetry, measurement 
or some non-perturbative method, and the remaining factors by perturbation theory up to 
0(a h (Q 2 )). We now wish to use the value of (f\02 k \^)\i} defined in this way to make a 
prediction for another physical process, 1Zfi(Q 2 ) sayQ, up to and including corrections of 

0(1/0"): 

n ft (Q 2 ) = d^q 2 /^ 2 ) ( f 1 1 i ) + D2iQ ^ 2) ( f I om !<> + •••. (3) 



4 Although this may present additional difficulties in practice. 

5 For simplicity we assume here that the scale Q 2 is the same in both processes; see below for the general 
case. 
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Returning to the example of R e + e - (Q 2 ) considered before, this corresponds to extracting the 
value of the gluon condensate from the cross-section at charmonium energies (the process 
V), and using it together with parton model perturbative calculations to make predictions 
for other processes. If we know D\ and D 2 at order k, then using eq.(^|) as the definition for 
{f\0^\ix)\i) we may write 



n fl (Q 2 ) = D?(f\0 1 \ i ) + ^ r (f\Oi k) \ i ) + 



D\^> - C\ 



(>k) n (>k)D 2 



(k) 



a 



(k) 



'a k + v 



(f\Oi\i)+o\-^y (4) 



where we have suppressed the renormalization scale (/x) and the arguments of the coefficient 
functions and of a s . The superscript (> k) denotes that the perturbation series for the 
coefficient function starts from the 0(a k+l ) term. Each of the terms on the r.h.s. of eq.(^) 
is free of renormalon ambiguities: the coefficients D± and D% because they correspond 
to series that have been truncated at a finite order in a s ; the matrix element (f\0^\n)\i) 
because it has been defined directly from an experimental measurement using eq.(0); and 
the coefficient in the third term because the renormalon ambiguity cancels in the difference 
of the two series up to a precision of 0(1/Q n+P ) (the perturbation series for the coefficient 
function D\ also has a renormalon ambiguity of 0(1/Q n ), which is cancelled by that in the 
series for C\(D2 /Czf). The low order terms in these two series are in general very different, 
but the divergent behaviour at high orders is controlled by the same renormalon singularity 
and hence is the same. 

The traditional procedure is to assume implicitly that the third term on the r.h.s. of eq.(^) is 
small, and to neglect it. For the small values of k for which perturbative results are generally 
available, there is in general no guarantee that the third term cannot give a contribution 
which is comparable to, or even larger than, that of the condensate itself. Moreover, the 
relative size of the two contributions will depend on the order k, the scale /i, and the process 
'PfiiQ 2 ) used to define the condensate, as well as on the process 1Zfi{Q 2 ) for which we want 
to make the prediction. 

We have seen that in order for the prediction for lZfi(Q 2 ) to be accurate up to and including 
terms of 0(1/ Q n ), the perturbative series for the coefficient functions must be known up to 
a sufficiently high order k, so as to make the third term on the r.h.s. of eq.(^) negligible. 
Formally, each term in these series is exponentially larger than the power corrections which 
are being calculated, and it may take a large number of terms before the required precision 
is achieved. In sections |||-|5| we investigate the accuracy of the theoretical prediction as a 
function of the order in several simple cases. 

In the above discussion we have assumed for simplicity that the two processes Vfi(Q 2 ) and 
lZfi(Q 2 ) occur at the same large momentum scale Q 2 []. This clearly is not necessary. For 

6 We also assume that the coefficient functions for the two processes are known at the same order k. 
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example one might wish to predict the behaviour of the moments of deep inelastic structure 
functions with the photon's momentum Q 2 , using an experimental measurement at a single 
value of Q 2 to determine the matrix element of O2. Similar cancellations of renormalon 
ambiguities occur also in these cases, because the Q 2 dependence of the coefficient functions 
is given by perturbation theory. 



2.2 Non-Perturbative Computation of the Matrix Elements of O2 

In this subsection, the procedure needed to compute the matrix elements of the higher-twist 
operator O2 (non-perturbatively) is described. This requires the operators O2 to be defined 
using a hard (dimensionful) ultraviolet cut-off A. For example, in lattice simulations it is 
natural to use bare operators, defined by the lattice action and with A = a -1 , where a is the 
lattice spacing. We therefore present the corresponding discussion in terms of bare operators, 
so that the renormalization scale \i is replaced by the cut-off A. We also have in mind that 
A 2 <C Q 2 , otherwise the use of the operator product expansion would not be necessary, one 
could just compute Vfi(Q 2 ) directly []. With a hard cut-off the operator O2 mixes with Oi, 
with mixing coefficients which diverge as A™. Thus, for the process Vfi(Q 2 ), we have 

VdQ 2 ) = c[ k \Q 2 /A 2 ) ( f I 1 (A) I i ) + g 2 fc) (QVA 2 ) ( ; ! 0a(A) ! ih (5) 

where the perturbation series, up to order k, for the coefficient function C\ takes the form 

C? \Q 2 /A 2 ) = cf\Q 2 /K 2 ) + cf } (Q 2 /A 2 ) (-)" (6) 

and the elements of the series C\ and Z\ diverge at most as powers of log(Q 2 /A 2 ). The term 
proportional to Z\ arises as a result of the mixing of O2 with 0\ |J. The matrix elements 
of 2 , computed with the hard cut-off, such as the lattice spacing, are well defined and 
unambiguous. Thus the same must be true for Ci, and the series in eq.© is indeed free of 
ambiguities of 0(1/Q n ). However this arises as a cancellation of the renormalon singularities 
in the two series cj[ an d 5i as k becomes large. Each of these two series has a renormalon 
ambiguity of 0(1/Q n ). The cancellation of this ambiguity occurs between contributions 
which, in each order of perturbation theory, are of different order in 1/Q. Again the low 
order terms in the two series are very different from each other, but the high order behaviour 
is governed by the same renormalon singularity and is the same. In order to predict Vfi(Q 2 ), 
we have to control the two series in eq.@ up to a precision of better than 0(AQ CD /Q n ). 

A related practical problem is the cancellation of the power divergences of 0(A n ) in the co- 
efficient function C\ with those in the matrix element ( / | 02(A) C\ has to be computed 

Although, in some cases this may not be possible in Euclidean space. 
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to sufficiently high order so that in eq.(|5|) this cancellation also occurs with a precision of 
better than 0(A% CD /Q n ). 

In addition to the ultraviolet cut-off necessary to regularize the theory, it is useful in some 
processes to introduce a physical cut-off. This is the case, for example, in the zero-recoil 



sum rules || |10fl , where the cut-off, A, is defined to suppress the contribution of states with 



excitation energies greater than A Q . There is a parallel between the above discussion with 
the hard cut-off A, introduced for the regularization, and that with the physical cut-off A. 
In the latter case, when, for example, the matrix element ( / | 02(A) | i ) is the kinetic energy 
Ai (which is of 0(Aq CD ) ), the corresponding power correction to the coefficient function (i.e. 

the term proportional to /Q 2 in eq.(|J)) contains perturbative corrections of 0(a s (A) A 2 ). 
The difficulty in achieving an accurate determination or bound with a physical cut-off follows 
the same discussion as for the hard ultraviolet cut-off. Only by arriving at an order such 
that a"(A) < 0(Ai/A 2 ), can one obtain significant results. 



2.3 Summary 

Before concluding this section, we briefly summarize the main points of the above discussion. 
In order to evaluate the power corrections to hard scattering processes, it is necessary to 
determine the matrix elements of the higher-twist operators, such as O2 in the above exam- 
ples. This can be done by comparing a theoretical prediction which depends on a matrix 
element of O2 to the experimenta data, or by computing the matrix element of O2 using some 
non-perturbative method. We have argued that, in the predictions for physical quantities, 
renormalon ambiguities, which are of the same order as the power corrections being evalu- 
ated, cancel in the combinations of coefficient functions given in eqs.@ and @. Either of 
these procedures reduces the "intrinsic" error of the calculation from 0(1/Q n ) to 0(1/Q n+P ), 
where by intrinsic we mean the minimum error achievable in principle. In order to reach the 
required precision, however, the series need to be evaluated to a sufficiently high order. An 
indication of whether the order is sufficiently high is given by the "common-sense" criterion 
that the last known term of the perturbative series of the leading coefficient functions should 
be significantly smaller than the power corrections. This point will be further illustrated in 
the next section, using a toy example, and in sections |] and [5], where the presence and can- 
cellation of renormalon ambiguities in quantities depending on the mass of the heavy quark 
and on the gluon condensate will be discussed. 



8 Note that A 3> Aqcd in order to be able to use perturbation theory in the parton sector; on the other 
hand we want A to be as small as possible in order to suppress the contribution of excited states and to 
reduce the uncertainty clue to the perturbative contributions which arc quadratic in A (sec below). 
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3 A Toy Model 



In this section we study a simple example, which contains many of the general features 
expected in operator product expansions, including the next-to-leading twist contributions, 
and the corresponding renormalon singularity. In this example, the "physical" quantity V 
is defined by 

V(x,e) = 10 3 due- u/x f(u), (7) 
Jo 

where 

1 f (l-2«) 1 

f(u) = ; ; (8) 

M ; l-2u l-2ur(l + 2u) ' K ' 

f(u) is the Borel transform of V(x, e), and e is a function of e and x defined in eq. fllOD below. 
The factor of 10 3 in eq.(|7[) is introduced for convenience. In eq.(|8]) the first term, expanded 
in powers of u and integrated over u as in eq.(|7|), generates the perturbation series of the 
leading coefficient function in x 



/ due~ u/x -> / due- u/x Y(2u) j 

Jo 1 — 2u Jo 



3=0 

x + 2x 2 + --- + ^-{2xf + ■■■ . (9) 



The second term in eq.(|8|) corresponds to the matrix element of the higher- twist operator, 
and we shall refer to its contribution to V as the "condensate" contribution. Its strength is 
governed by the parameter e, defined by 

e = ee- 1/2x } (10) 

which shows that e is a term of order Aq CD /Q when the "coupling" x ~ 1/ ln(Q 2 /AQ CD ) []. 
In physical cases, e is determined by the non-perturbative dynamics, so that V is only a 
function of x. In our toy model, however, we will treat e as a free parameter. Its natural 
value is of 0(1). The factor of 1/T(1 + 2u) has been introduced in order to make the integral 
in eq. (0) converge at large values of u (for small values of e) or to improve the convergence. 
Other choices of such damping factors at large u would have been equally good for our 
purposes. 

Both terms in eq.(H) exhibit a renormalon singularity at u — 1/2; this however cancels in 
f(u), in a way which is analogous to the cancellation that is expected to occur in physical 
cases. Because of this singularity, the coefficients of the perturbation theory grow like a 
factorial. As a consequence, even for small values of x, the contribution of high orders 

9 Similar models can be constructed to mimic the cases in which the power corrections are suppressed as 
(Aqcd/Q)", with n > 1. 
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e = 100 


e = 10 


6=1 


e = 0.1 


V{x = 0.07, e) 


72.23 


82.33 


84.92 


199.57 


k 


tk 


Sk 


"Value of Condensate" - 


C k (x,e) 


1 


70 


70 


2.23 


12.33 


14.92 


129.57 


2 


9.8 


79.8 


- 7.57 


2.53 


5.12 


119.77 


3 


2.74 


82.54 


- 10.32 


-0.21 


2.38 


117.02 


4 


1.15 


83.70 


- 11.47 


-1.36 


1.23 


115.87 


5 


0.65 


84.34 


- 12.12 


-2.01 


0.58 


115.23 


6 


0.45 


84.79 


- 12.57 


-2.46 


0.13 


114.77 


7 


0.38 


85.17 


- 12.95 


-2.84 


-0.25 


114.39 


8 


0.37 


85.55 


- 13.32 


-3.21 


-0.62 


114.02 



Table 1: Perturbative and non-perturbative contributions to V, for x = 0.07. The full result 
for V is given in the second row. and Sk are the k-th term and the sum of the first k terms 
of the perturbation series @. The condensate contributions Ck are defined in eq. (|TTJ) . 



diverges. We denote by k m i n the order in x for which the magnitude of the k-th term of the 
series is the smallest one; k min depends on x only. 

In order to mimic the procedure one is forced to adopt in realistic cases, we proceed as 
follows: 

i) For any choice of x and e, the value of V(x,e) is obtained exactly from eq.(|7|) by 
numerical integration, and represents the "experimental" result (with no error). 

ii) We assume that only k terms, with k < k min , of the perturbation series (H) are known. 

iii) We define the "condensate" Ck as 

C k (x,e)=V(x,e)-lj2r(j)(2xy , (11) 

where by writing explicitly the arguments of Ck we recall that e is being treated as 
a free parameter. The subscript A; is a reminder that the condensate was defined by 
subtracting k terms of the perturbation series. 

iv) It is then envisaged that the values of the condensate obtained in this way are used to 



make predictions for another process, 1Z say, as explained in subsection 2.1 



The stability of the value of Ck with k, and the comparison of its value to the k-th term of the 
perturbation series, monitors the precision that can be reached in predictions for physical 
quantities containing corrections of 0(Aqcd/Q)- As an illustration we present in table |l] 
the numerical results for x = 0.07 (for which k m i n = 8) and for e = 100, 10, 1 and 0.1 . The 
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"physical" value of V[x = 0.07, e) for each of the values of e is given in the second row of 
the table, tk is the k-th term of the perturbation series @ and Sk is the sum of the first k 
terms, and they are tabulated up to k = k min = 8. The remaining entries in the table are 
the values of the "condensates", Ck, which one would deduce at each order of perturbation 
theory. We now comment on these, distinguishing between large, intermediate and small 
values of e in turn: 

i) Large value of e (e = 100): 

For e = 100, the values of the condensate (Ck) stabilize quickly, so that already for 
k = 3, C3 is significantly larger than the corresponding term in perturbation theory 
(1(3) and the sum of the higher order terms from t 4 to t^ min . If this is the case for both 
V and 1Z (in the notation of section |2.1|) then the precision of the prediction for 1Z is 



clearly improved by including the contribution from the condensate, provided that at 
least three terms in the perturbation series have been calculated. 

ii) Intermediate value of e (e =10 and e = 1): 

As e is decreased, one has to calculate more terms of the perturbation theory before 
the values of Ck stabilize. For e = 10, the improvement in including the condensate 
contribution is, at best, marginal, even if five or six terms of the perturbation theory 
have been computed. For e = 1, V is estimated accurately by perturbation theory, 
and the condensate contribution is too small to be determined. In many practical 
situations it may be sufficient to know that the condensate is smaller than some value. 
The uncertainty in the perturbation series is independent of the value of e. 

iii) Small value of e (e = 0.1): 

For small values of e the situation becomes very unstable, and the contribution of per- 
turbation theory, even if one includes all the terms up to tk min , is a poor approximation 
to V. The reason can be understood by considering the series in x generated by each 
of the two terms in eq.(§). The coefficients of each of these two series grow like the 
factorial of k, but this factorial growth is cancelled in their sum. For very small values 
of e, however, this cancellation will only begin to take effect at very high orders, in 
particular at values of k such that k 3> k min , at which the perturbation series is 
already diverging rapidly with k. It would be fascinating to find a realistic physical 
example corresponding to this case. 

The conclusions which we draw from this simple example are as follows. Imagine that 
V(x, e) has been measured, and that k terms of the perturbation series have been calculated. 
A reasonable estimate of the perturbative contribution to V would then be Sk ± tk- If 
Ck tk, then it makes sense to call Ck the condensate contribution, and to use it in 
predictions for other processes for which the perturbation series has similar properties and 
has been calculated to the same precision. The challenge in phenomenological applications 
is to demonstrate that this is the case. Otherwise the condensate is of the same order as the 
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uncertainty in the perturbation series (or smaller) and hence its value would depend on the 
process from which it is extracted, tk serves as an estimate of the uncertainty in the value 
of the condensate. 

There is one further important point which we wish to stress, i.e. the universality of the 
higher order corrections. Consider a set of processes V, TZ ■ ■ ■ for which the leading power 
correction is given by the same matrix element of a given higher-dimensional operator. The 
high order behaviour of all the leading coefficient functions is then dominated by the same 
infra-red renormalon, and is hence universal. This leads to the possibility that the error in 
the prediction for process 7Z, obtained by using the condensate Ck determined in process V, 
may be smaller than the estimated uncertainty in the value of the condensate itself. In 
order to see this, note that the definition of the condensate in eq.([TT|) implies that 

tfs+i + ^fe+i ~ — , (12) 

where the superscript on t^ +1 implies that this is the (&+l)-th term in the perturbation series 
for the process V. The difference between the (k + l)-th order and fc-th order predictions 
for TZ is 

tf +1 + C k+l -C k + (*^p) , (13) 

where the last term represents perturbative, renormalon-free corrections to the power sup- 
pressed term. If the perturbation series for the two processes are both dominated by the 
same renormalon, and k is sufficiently large, then it may be that ~ t£ +1 , and that the 
(k + l)-th contribution to the prediction for TZ is smaller than the difference in the conden- 
sates, Ck+i — Ck- In practice, however, often only one or two terms of the perturbation series 
are known and it is unclear to what extent this property of universality will be useful in 
phenomenological applications. 



4 The Pole Mass of a Heavy Quark 

In this section the discussion of section || is applied to the computation of the heavy-quark 
mass in the HQET. Of course in practice we do not know the perturbation series for the 
coefficient functions to sufficiently high order to be able to study the numerical effects of 
the cancellation of renormalon ambiguities directly. For this reason we present a calculation 
performed in the limit of a large number of light quark flavours (Nf), or more precisely we 
perform the perturbative calculations keeping only the term with the highest power of (3q at 
each order: (3q = 11 — 2/3Nf is the lowest order coefficient in the /3-function. 

It is possible to compute the renormalized mass of a heavy quark, defined at a large renor- 
malization scale, \x 3> Aqcd, from the matrix elements of the HQET obtained with some 
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non-perturbative method, such as lattice simulations |TT1. As an example we consider m, 



the mass defined in the MS scheme at a renormalization scale /i = m: 

m = m m (m m ) . (14) 

The computation requires the expansion of the propagator of the heavy quark in QCD (with 
the mass, wave-function renormalization and coupling constant defined in the MS scheme 
say), in terms of matrix elements of operators in the HQET. We will see that the coefficient 
function of the leading operator, C\ of eq.(|l|), is indeed a difference of two series (as in eq.(||)), 
each of which has a renormalon ambiguity of 0(Aqcd) ; the ambiguity cancelling in the 
difference. In the large Nf limit the terms of the series can be calculated to arbitrarily high 
orders, and the cancellation of the ambiguity then observed. This is done in subsection 



in subsection [4.1| we start with a brief review of how one computes m using non-perturbative 



methods such as lattice simulations. 



4.1 Evaluation of m in Lattice Simulations 

In this subsection the procedure needed to evaluate m from simulations in the HQET up to, 
and including, terms of 0(Aq CD ), but neglecting terms of 0(AQ CD /m) is briefly reviewed 
HTT| 0. The non-perturbative quantity which is computed directly in lattice simulations, and 
which is required for the determination of m, is the bare binding energy where the label 
H denotes the hadron containing the heavy quark. Sh is obtained from the time dependence 
of the correlation function of two interpolating operators (Jh) for the hadron H: 

^2(0\J H {x,t) 4(0,0) | 0) = Z e~ £Ht , (15) 

X 

where t is sufficiently large for the correlation function to be dominated by the lightest 
particle created by J H , which is assumed to be H. 

The relation between Sh and m can be obtained by matching the heavy quark propagator in 
QCD with operator matrix elements evaluated in the HQET fLlf . Using the MS renormaliza- 



tion scheme at a scale /i for the mass, wave function and coupling constant renormalization, 
the inverse propagator in QCD {S~ l ) is of the form: 



n=0 



+v-kjr d n (v ■ k/fi, m{tiln) + 0(A 2 QCD /m) (16) 



n=0 



10 In ref.|Tl| the generalization of this discussion to include terms of 0(Aq CD /to) is also presented. 
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where the momentum of the heavy quark is tjiqv + k v is the four velocity of the heavy 
quark and Sp is defined by 

l±i Sp = l±l s l±i. (17) 

2 2 2 K J 

The first series on the right-hand side of eq.([T6|) is just the perturbative expansion of the 
pole mass of the heavy quark in terms of m(/x): 

"Vie = m(n) [^^J c n (m(n)/n) . (18) 
The Borel transform of the series in eq.flTgp has renormalon singularities, which is a ma- 



nifestation of the fact that the pole mass is not a physical quantity || ||. The ambiguity 
corresponding to the leading singularity is of 0(Aq CD ). 

Now consider perturbation theory in the HQET (defined by the action hv-D h), using the 
lattice spacing as the ultraviolet cut-off. Sp 1 can be expressed in terms of S~$, the propagator 
in the HQET, which can be interpreted as the matrix element of the operator O2 = hv-D h: 

S^iv-k) = m Q - (m po ie - 5m) 

+ C v . D (m(fi)a,m(fi)/fi)S~g(v ■ ka) + 0(A 2 QCD /m([i)) , (19) 



fo= ^ff«Mr^ ,20) 



where the series 

= 1£= v (<- 

Arc 

it — u \ J 

is just the perturbative expansion of the inverse of the quark propagator in the HQET at zero 
momentum. The leading renormalon singularity in m po i e is cancelled by the one in the series 
for 5m, so that the combination m po i c — 5m has no renormalon ambiguity of O(Aqcd) H- 
The two series in eqs. (|T8"D and are an example of the series in eq.(^). Here A = a -1 , Q 
is the mass of the heavy quark and n = 1, i.e. we are calculating the O(Aqcd) correction to 
the mass of the heavy quark, which is one power of m smaller than the leading term. The 
operator 0\ is hh, which is a conserved current in the HQET with matrix element equal to 
1, and so it does not appear explicitly in eq.(]l9"D. 

The mass m is obtained from the relation 

"Vic - Sm = M H - £ H , (21) 



where M# is the physical mass of the hadron H, and by inverting the relation in eq. (|18|) 
between the pole and MS masses. In eq . (plf) , the linear divergence present in Eh is cancelled 
by that in the series 5m, and the renormalon in m po i e is cancelled by that in 5m, as explained 
above. In the following subsection we study the numerical cancellation of the renormalon 
ambiguity between m po i e and 5m. Numerical results for m obtained in this way (but with 
the perturbative terms only computed to one- loop order) have been presented in |L1|, [L2]. 

11 The precise definition of toq here can be conveniently chosen later. 
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Figure 1: Diagrams which must be evaluated in order to study the heavy quark mass (and 
hence A) in the large-/?o limit. The double line represents the propagator of the heavy quark, 
and the bubbles represent light-quark loops. 




Figure 2: One-loop bubble graph contributing to the vacuum polarization of the gluon. This 
graph is the basic ingredient in the evaluation of the large-/?o contribution to the quark mass 
(and other quantities). The solid lines represent light-quark propagators. 



4.2 Cancellation of Renormalon Ambiguities in the Heavy-Quark 
Mass 

In this subsection we explicitly trace how the cancellation of the renormalon ambiguities 
occurs in the combination m po \ e — 8m. The calculation is performed in the large-/?o limit, in 
which only the terms containing the leading power of /3q are kept in each order of perturbation 
theory. As a further simplification, we perform this calculation in the large-/?o limit, with 
the Pauli-Villars (PV) cut-off A as the ultraviolet regulator in the effective theory (rather 
than the lattice spacing) 0. In other words, in analogy with lattice field theory, we imagine 
that we have computed Sh non-perturbatively in the PV theory for some hadron H, and 
now perform the calculation of the matching term m po i e — 5m in perturbation theory. 

In order to obtain the result in the lowest non-trivial order in the large-/?o approximation, 
it is sufficient to evaluate the set of diagrams in fig. [I], summing the contributions from an 
arbitrary number of light-quark loops. Consider a single light-quark loop insertion as in 
fig. |2|. We denote the expression from this diagram by L(p,m^), where m q is the mass of 
the light quark. L is given in terms of a divergent integral, which we regulate by taking 
L(p, 0) + L(p, 2A 2 ) — 2L(p, A 2 ) as the regulated expression for the diagram, as suggested in 
ref.[|]l|. Other choices are also possible. For the gluon propagator it is sufficient to take 

12 The general features of the cancellation of the renormalon ambiguities are qualitatively the same with 
any hard cut-off. For the purposes of illustration the PV cut-off is very convenient. 
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—ig^iXjq 2 — l/(q 2 — A 2 )) since any terms proportional to q^q v in the gluon propagator do 
not contribute to 5m. We then find 

where 

/•OO rjrr, 

K n = I (1 + x2)n+1 [5/3 - ln(2* 2 ) - F(x 2 /2) + 2F(x 2 )} n , (23) 
a s (A) is the bare coupling constant and the function 

F(x 2 ) = 6 f dyy{\ - y) ln(l + y(l - y)x 2 ) (24) 

can be readily evaluated. The position of the leading infra-red renormalon singularity is 
known, and by inverting the Borel transform using the saddle point method we find that the 
behaviour of K n at large n is given by: 

5/6 

K n - ^ n n\ (25) 

The numerical results obtained by using eq. (|23|) approximate their asymptotic values in (|25| ) 
to better than 5% already for n — 4. 

The relation between the pole mass and m in the large-/?o limit is given by |TJ| : 



"^poio = m 



(26) 



where ac s (m) is the MS coupling constant at the renormalization scale m. For large values 
of n, the coefficients d n behave as 

d n -> e 5/6 2 n n! (27) 

and again this asymptotic relation is well satisfied by the numerical results for relatively 
small values of n. 

To demonstrate the numerical cancellation of the leading renormalon singularity, we do the 
following. In the second column of table § we present the results for the pole mass obtained 
from eq. (]2"6j) , with the assumption that m = 4.5 GeV, and a s (m) = 0.2. The results are 
presented in successive orders of perturbation theory in the MS coupling constant at \i — m. 
As expected, the value of m po i e increases rapidly at high orders due to the presence of the 
renormalon singularity. The standard approach when dealing with an asymptotic series is 
to consider only the first few terms to estimate the result. Since the smallest term in the 
series is of 0(100 Me V) (for n = 3 — 5), i.e. of O(Aqcd) as expected, this can be viewed 
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Tfl i C-rPV 

" ''pole v 


777 i — A 777 rip\f 
pole \j 1 1 v v 


tree 


4.5 


4.5 





4.88 


4.62 


1 


5.12 


4.72 


2 


5.24 


4.74 


3 


5.34 


4.76 


4 


5.44 


4.77 


5 


5.58 


4.77 


6 


5.80 


4.78 


7 


6.21 


4.78 


8 


7.08 


4.79 



Table 2: m pole and m po i — 5m calculated up to 0(a™ +1 (m)) in perturbation theory, in the 
large-/? limit, using m = 4.5 GeV and a s (m) = 0.2. 



as the intrinsic uncertainty in the unphysical quantity m po i e . In the third column of table ^| 
we present the analogous results for m po \ e — 5m, also expanded in terms of the MS coupling 
constant at fi = fn. In eq. (^) we have taken A = 2 GeV, and have expanded the Pauli-Villars 
coupling a s (A) in terms of the MS coupling a s (m). The series of the difference m po i e — 5m 
stabilizes at lower orders because of the cancellation of the leading renormalon singularity. 
The remaining uncertainty is now much smaller (in general it would be of 0(AQ CD /m), but 



in the large-/?o limit the corresponding renormalon is absent [§], [15|] , so that the ambiguity is 
of 0(AQ CD /m 2 )). As a consequence the smallest term in the series, which sets the scale of 
the intrinsic uncertainty, is now less than 10 MeV. 

The unphysical parameter A = Mh — rn po i e is frequently used in phenomenological studies 
of B-physics. Results and bounds for A are presented (see for example |16| and the reviews 
|i~7| , [i~8| , and references therein). In order for this to make any sense, A must be defined 
precisely in terms of some physical quantity, and its value will depend on this physical 
quantity and on the order of perturbation theory used to extract A. As can be seen in 
table |2|, the values for A will change by several hundred MeV as the order of perturbation 
theory is increased. This also implies that when using the value of A measured in one process 
to make predictions for a second one, we should use the same order of perturbation theory 
in both processes and hope that the universality discussed at the end of section [| holds to a 
good approximation. 

All calculations of power corrections to hard scattering and decay processes will involve 
a cancellation of renormalon ambiguities similar to the one discussed in this section 0. Of 
course, one should remember that the results in table |^ were obtained using an approximation 
and can only be taken as being indicative. Nevertheless they highlight the difficulty of 



13 Unless there is a symmetry which prevents the mixing of the operators with different dimensions. 
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evaluating power corrections. If one's aim is to try to evaluate m up to uncertainties of 
0(Aq CD /m) ~ 25 MeV or so, many orders of perturbation theory would be required because 
the series for m po i e — 8m converges very slowly. In the lattice theory only the one-loop term 
is known (corresponding to n = in table |2]), and the continuum relation (|18|) is known 
up to two-loop order. This would suggest that values for the binding energy A which are 
used in phenomenological studies in heavy-quark physics are uncertain by an amount of 
order 100 MeV, due to our ignorance of the higher order perturbative terms. A similar 
comment applies to our computation of m in ref . [|TTJ where this uncertainty was probably 
underestimated, and to the result in |T2]] . 

It is likely that the case of A is a relatively good one, since the corrections are suppressed by 
just one power of the mass of the heavy quark, and, at least in the large-/?o approximation, 
the asymptotic behaviour ( P5"| ) and Q2"7| ) seems to set in at low orders of perturbation theory. 
In the next section we consider an important example where this is not the case, that of the 
evaluation of the gluon condensate. 



5 The Gluon Condensate 

In this section we study another important example, that of the contribution of the gluon 
condensate to physical quantities in general, and to the D-function in e + e~ annihilation in 
particular. The D-function is defined by D{Q 2 ) = —\/AQ 2 dH{Q 2 )/dQ 2 ) where II is obtained 
from the correlation function of two electromagnetic currents: 

% J d*xe^ x (0\T{J^(x)J u (0)}\0) = (q,q u - g^q 2 )U(-q 2 ) (28) 

and Q 2 = —q 2 . The gluon condensate, (a s G 2 /n) , is the vacuum expectation value of 
an operator of dimension 4, and hence its contribution is of 0(Aq CD /Q 4 ) relative to the 
perturbative terms. As we have tried to stress throughout this paper, the gluon condensate 
itself is not a physical quantity, as it contains a renormalon ambiguity of 0(Aq CD ). This 
ambiguity is cancelled by that in the perturbation series for the D-function (or, in general, 
by that in the leading coefficient function for the process being studied). The suppression 
by four powers of Q implies that the cancellations are very large and leads to enormous 
difficulties in the quantitative evaluation of the power corrections. 



In subsection ^2| below we study the high order behaviour of the perturbative series for 
the D-function in the large (3q limit. We argue that even at low values of Q 2 , where the 
relative contribution of the power corrections is significant, the uncertainty and ambiguity 
in the perturbation series can be comparable to the contribution normally ascribed to the 
condensate. We start, however, by a discussion of the computation of the gluon condensate in 
lattice simulations. We demonstrate that the extremely large numerical cancellations which 
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arise in the subtraction of the quartic power divergence (i.e. of the terms which diverge as 
a~ 4 ), and the presence of the corresponding renormalon singularity, make the quantitative 
evaluation of the leading power corrections (i.e. the 0(1/Q 4 ) corrections) prohibitatively 
difficult. 



5.1 Evaluation of the Gluon Condensate in Lattice Simulations 

A natural definition of the condensate in lattice QCD is given in terms of the expectation 
value of the plaquette variable P^ v (// and v define the plane containing the plaquette, but 
the expectation value is, of course, independent of the choice of plane): 

P = (l- -TrP, u ) = -a 4 ( —G 2 ) latt , (29) 
6 Sb it 

where a is the lattice spacing and the subscript stands for "lattice". The variable P is 
measured very precisely in lattice simulations for the Wilson action, at all standard values 
of the lattice spacing (as well as for some other lattice discretizations of QCD). In lattice 
QCD, as with any regularization using a hard cut-off, P is not zero in perturbation theory, 
but is given by an expansion of the form 

^ = £^> (30) 

n=l P 



where (3 = 6/<7g(a) and go(a) is the bare lattice coupling constant. The series in eq.(|30D arises 



as a result of the mixing of the G 2 with the identity operator. The first 8 (!) coefficients q 
have been obtained numerically using Langevin techniques ||19 |. We would like to use the 



computed value of P and perturbative matching to calculate the 0(Aq CD /Q 4 ) corrections 
to some physical process. For example for the D-function the relation is: 



D(Q 2 ) = D£ - + (l + I +~')P, (31) 

where 



and 



D cont = 1 + ^M + 1M0 (^M\ +... (32) 




Dp°°* is the perturbative series for the D-function and the superscript stands for "contin- 
uum" . The series D l ^ t arises from the matching of the D-function with the lattice operator 
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in eg. fl29|) , and the superscript stands for "lattice". In the notation of subsection [2.2| , and 
eqs. (HD and (||) in particular, 

D c ™± = Cl (A 2 /Q 2 ), (34) 



and 



^pert = -Ci(Q M ) I g ) • (3.-)) 



with A = a 1 and n = 4. In this case the operator 0\ is the identity operator. The leading 

p°rt is cancelled by the renormalon in D^ t , 



infra-red renormalon singularity in the series is cancelled by the renormalon in _D lat 



as explained in subsection |2T2| . 

As an example of how serious the cancellations are, and how very difficult it is to obtain 
results with the required precision we take Q = m T {m T is the mass of the r-lepton), 
(3 = 5.7, which corresponds to an inverse lattice spacing of about 1.15 GeV, and af s (m T ) 
= 0.32. We are forced to choose fairly large values of Q 2 , for which the contribution of the 
condensate is expected to be small, since in order to make use of the lattice results we require 
Aqcd < a -1 < Q. We consider (3 = 5.7 (a -1 ~ 1.15 GeV) to be about the smallest value 
of P (and inverse lattice spacing) for which one may reasonably expect that lattice artefacts 
will not invalidate the interpretation of the results. For these values of the parameters, the 
three components in eq.fl3~H) begin like: 

L>cont = 1 + 0.102 + 0.017 + ■■■ (36) 
l£tt = 5 630 _ 6.259 + - •• (37) 

24 / 7 ry m (Q) \ 

' 1 + - 8 W) + • • • }P = (1 + 0.119 + •••) 1.894 . (38) 



a 4 Q 4 V 6 7T 



The numbers in eq. fl37|) have been obtained after rewriting the series in terms of af s (Q). 
We see that the terms in eqs.(|3"T|) and ( |3"ED are huge compared to the contribution one would 
normally ascribe to the gluon condensate of about 1% or so (taking ( a s G 2 / ir) ~ 0.018 GeV 4 ) . 
It would clearly be enormously difficult to quantify these power corrections accurately. This 
would require the perturbation series -Dpert — -Dpert; an d also the relation between the lattice 
and MS coupling constants, to be known to extremely high orders. 

Although the numerical results presented above were obtained using the plaquette variable 
to define the gluon condensate on the lattice as in eq. (|29|) , the general discussion applies to 
any choice of operator. The problems arise from the presence of quartic power divergences 
and renormalon ambiguities, which are general features of lattice attempts to evaluate the 
power corrections associated with the gluon condensate. 

In ref. |^0[ the author has attempted to define the gluon condensate from P — J2 c n/P n , 
where P has been measured numerically, and the first 8 coefficients c« are known [113. He 



uses different resummation techniques for the perturbative terms and finds results for the 
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condensate which depend significantly on the method of summation, and are always at least 
five times larger than those used in phenomenological applications. The point that we are 
trying to stress in this paper is that such an analysis is theoretically inconsistent. The series 
J2 c n//3 n has a renormalon ambiguity of 0(Aq CD ) (which, as always, is of the same order as 
the effect one is trying to evaluate); thus subtracting the perturbative series J2 c n/P n from 
the computed values of P leaves an intrinsic arbitrariness of this order. In order to eliminate 
this arbitrariness we have to include the measured value of P in the prediction for a physical 
quantity using the matching procedure described above. In this way, for the D-function 
(which is a typical example) the renormalon ambiguity cancels between the first two terms 
on the right-hand side of eq. (|31|). 



5.2 The Uncertainty in the Perturbation Series for the D-function 

Although the discussion in the previous subsection was concerned specifically with the eval- 
uation of the power corrections to the D-function using lattice simulations, we believe that 
it is also very difficult to control the corresponding calculations in phenomenological studies 
using continuum regularizations. Consider the perturbation series for the D-function using 
the MS coupling constant: 

£cont = 1 + aJ&l + L639g fajgy\ + g 371Q1 faj^)\ + . . . . (39) 



The coefficient of (a s (Q 2 ) /n) 4 has been estimated to be about 27.5 PJ], based on calculations 



using the principle of minimal sensitivity p2] and the effective charge approach |23]. For 



values of Q 2 such that a s {Q 2 ) < 1/2 say, the series appears to be reasonably well behaved. 
However, in order to gain some insight into the effects of the infra-red renormalon in the series 
-Dp°rt we have to go beyond the order for which the coefficients are known. For this reason 
we study this series in the large-/?o limit, for which the coefficients have been determined in 
refs. |24|, |25|. Writing 

fl -, + 'fJWfl)*, (40 ) 

PO n= l V 7T / 

the coefficients readily be obtained from the Borel transform of the series, 

~ (Au) n - 1 32e~ Cu ~ (-l) k k 

h r(n) ~ 3(2 - U ) jfr 2 [fc2_(i_ u) 2]2 ' I ) 

where in the MS scheme C = —5/3. The large order behaviour of the series is dominated 
by the singularity closest to the origin, which in this case is a double (ultraviolet) renormalon 
pole at u — — 1. The high order terms generated by this pole diverge like a factorial of the 
order, but with alternating signs, corresponding to a behaviour which is Borel-summable. 
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n 











l 


1 


1 


1 


1.049 


0.049 


4.461 


2 


1.089 


0.040 


0.799 


3 


1.125 


0.036 


0.286 


4 


1.161 


0.036 


0.154 


5 


1.201 


0.039 


0.110 


6 


1.248 


0.047 


0.099 


7 


1.310 


0.062 


0.106 


8 


1.399 


0.090 


0.133 


15 


23.511 


13.016 


13.536 



Table 3: Values of s n , t n and t ni asym P for a value of Q 2 such that a s (Q 2 ) = 1/2 (Q ~ 
0.75 GeV). 



Thus although the presence of ultraviolet renormalons may add further practical difficulties 
to the evaluation of the power corrections of 0((Aq CD /Q) 4 ), we will not consider them 
further here. Specifically, we subtract the contributions of the ultraviolet poles at u = — 1 
and u = —2, which appear in the terms with k = 2 and 3 respectively in eq. (41), by 
considering the behaviour of the coefficients n' n obtained from 

~ (AuY 1 - 1 , 32e~ CM f 90 - 39m + 5m 2 224 - 72m + 7m 2 ~ (-l) fc k 1 
r(n) Kn ~ 3(2 -u)\ 144(3 - u) 2 576(4 - u) 2 + ^ 4 [k 2 - (1 - u) 2 } 2 J ' 

(42) 

The residues of the infra-red renormalons (at u — 2 and above) are the same in eqs. ( flip and 
fl4"2]). The large order behaviour of the coefficients K* n is given by 

, e 10/3 

K n * K n,asymp gn— 1 ^ ' (^) 

asn-> oo. 

Consider the perturbation series generated by the coefficients n' n , -D™^* = J2 n =o where 
t = 1 and 

Po V ^ / 

for n > 1. We denote by s n the sum of the series up to ra-th order, s n = X^ =0 ifc. It is 
also convenient to define t niasymp as in eq.(|44"D, but with k' u replaced by the asymptotic form 
k' 

n„asymp - 

In order to illustrate the difficulties of evaluating the perturbation series with sufficient 
precision to make the inclusion of the corrections of order (Aqcd/Q) 4 meaningful, we present 
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an example. In table we give the values of s n , t n and t n ,asym P obtained for a value of 
Q 2 such that the MS coupling constant a s {Q 2 ) = 1/2 (Q ~ 0.75 GeV). We choose a small 
value of Q 2 so as to enhance the contribution of the power corrections (it would be even 
more difficult to quantify the power corrections at higher values of Q 2 ). We now make some 
comments on these results: 

i) The smallest term in the series {t n } occurs for n = 4, = 0.036. One might therefore 
be tempted to take the four-loop result, S4 = 1.161, as the best estimate for the sum 
of the perturbation theory, and £4 ~ 4% as the estimate of the renormalon ambiguity. 
One might also expect that the inclusion of the gluon condensate will eliminate this 
uncertainty (up to a precision of order (Aq CD /Q) 6 ). Moreover since with standard 
phenomenological values of the condensate ((a s G 2 /n) ~ 0.018 GeV 4 ) it is expected 
that its contribution to the D-function at such low values of Q 2 is about 40%, the 
ambiguity of about 4% can be considered negligible. 

Such an interpretation is wrong, however. The perturbation theory for n ~ 4 is not yet 
dominated by the leading infra-red renormalon (as can be seen, for example, from the 
fact that t 4 is very different from t^asymp, or by evaluating the contribution from the 
next-to leading renormalon at u = 3). The smallest term in the series {tn,asymp} occurs 
at n = 6 and is about 10%. Even for n ~ 6, however, the perturbation series {t n } is 
not well approximated by {t n ,asymp}- Thus it is not easy to estimate the uncertainty in 
the evaluation of the perturbation series, other than to say that it is certainly greater 
than 10%. 

ii) Since the perturbation series approaches its asymptotic value very slowly, there is 
no reason why the low order terms should be approximately universal. For exam- 
ple, the large order contribution to K,' n from the infra-red renormalon at u = 3 is 
— 4/27e 5 n!/12 n_1 and for n < 4 is greater than or comparable to that of the leading 
infra-red renormalon at u = 2 (which is e 10 / 3 (n — l)!/8 n ~ 1 , see eq.([13|) ). The renor- 
malon at u = 3 corresponds to operators of dimension 6, whose contribution relative 
to the gluon condensate depends on the process, and hence the low order terms of the 
perturbation series are not universal 0. 

iii) The discussion in i) and ii) was based on the assumption that it is possible to calculate 
many orders of perturbation theory and to study the extent to which the asymptotic 
behaviour has been reached. Of course, in practice, usually only one or two terms 
of the perturbation theory are known, which adds substantially to the uncertainty. 
For example the values of the contribution from the "condensate" , which one would 
obtain by subtracting either the one-loop result (si) or the six-loop one (sq) from the 

14 A related question is whether at such low scales, the non-perturbative contibutions of the operators 
corresponding to u — 3 and above do not become as large as that of the gluon condensate. 
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measured value of D(Q 2 ), would differ by about 0.2, i.e. by 20% of the D-function 
itself. 

iv) It may be the case that for some processes, and at small values of Q 2 in particular, 
the "condensate" contribution is much larger than all the combined uncertainties. In 
the language of subsection ^7T] this would be necessary both for the physical quantity 
V being used to determine the condensate (e.g. some correlation function used in the 
study of the spectrum of charmonium) and for the quantity TZ for which the prediction 
is being made (e.g. decay constants or semileptonic form factors of heavy mesons). 
This example demonstrates, however, that to be confident that the uncertainties are 
indeed sufficiently small will require considerable effort. 

In this simple example the uncertainty and ambiguity in the perturbation series for D(Q 2 ) 
is at least a significant fraction (50-100% ?) of the expected size of the leading power 
correction. This is in spite of the fact that Q 2 was chosen to be small in an attempt to 
minimize the relative size of the ambiguity. Although one can change the details of the 
discussion by using different values of Q 2 , different renormalized coupling constants as the 
expansion parameters, or different physical processes, it is our contention that the difficulties 
discussed above are general and cannot be easily overcome. 



6 Conclusions 

In this paper we have studied several examples in order to understand whether it is possible 
to compute power corrections to hard scattering and decay processes to a sufficient level 
of precision. By this we mean that the theoretical uncertainties must be smaller than the 
power corrections themselves. In all of the examples considered, we have found that this is 
not possible unless we are able to control perturbation theory at higher orders than those 
available at present. Since the arguments discussed in this paper are general, and not specific 
to the examples used, we believe that even the leading power corrections are currently not 
well determined, and there is little hope to compute higher order power corrections 0. We 
hope that these disappointing conclusions and provocative comments will help to stimulate 
further debate and a systematic investigation of this central question of particle physics 
phenomenology. 

15 For exceptions to the general discussion see section 0. 
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